Global arthropod beta-diversity is spatially and temporally structured by latitude.

Global biodiversity gradients are generally expected to reflect greater species replacement closer to the equator. However, empirical validation of global biodiversity gradients largely relies on vertebrates, plants, and other less diverse taxa. Here we assess the temporal and spatial dynamics of global arthropod biodiversity dynamics using a beta-diversity framework. Sampling includes 129 sampling sites whereby malaise traps are deployed to monitor temporal changes in arthropod communities. Overall, we encountered more than 150,000 unique barcode index numbers (BINs) (i.e. species proxies). We assess between site differences in community diversity using beta-diversity and the partitioned components of species replacement and richness difference. Global total beta-diversity (dissimilarity) increases with decreasing latitude, greater spatial distance and greater temporal distance. Species replacement and richness difference patterns vary across biogeographic regions. Our findings support long-standing, general expectations of global biodiversity patterns. However, we also show that the underlying processes driving patterns may be regionally linked.

Towards a toolkit for global insect biodiversity monitoring.

Insects are the most diverse group of animals on Earth, yet our knowledge of their diversity, ecology and population trends remains abysmally poor. Four major technological approaches are coming to fruition for use in insect monitoring and ecological research-molecular methods, computer vision, autonomous acoustic monitoring and radar-based remote sensing-each of which has seen major advances over the past years. Together, they have the potential to revolutionize insect ecology, and to make all-taxa, fine-grained insect monitoring feasible across the globe. So far, advances within and among technologies have largely taken place in isolation, and parallel efforts among projects have led to redundancy and a methodological sprawl; yet, given the commonalities in their goals and approaches, increased collaboration among projects and integration across technologies could provide unprecedented improvements in taxonomic and spatio-temporal resolution and coverage. This theme issue showcases recent developments and state-of-the-art applications of these technologies, and outlines the way forward regarding data processing, cost-effectiveness, meaningful trend analysis, technological integration and open data requirements. Together, these papers set the stage for the future of automated insect monitoring. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Decline in Honeybees and Its Consequences for Beekeepers and Crop Pollination in Western Nepal.

In understudied regions of the world, beekeeper records can provide valuable insights into changes in pollinator population trends. We conducted a questionnaire survey of 116 beekeepers in a mountainous area of Western Nepal, where the native honeybee Apis cerana cerana is kept as a managed bee. We complemented the survey with field data on insect-crop visitation, a household income survey, and an interview with a local lead beekeeper. In total, 76% of beekeepers reported declines in honeybees, while 86% and 78% reported declines in honey yield and number of beehives, respectively. Honey yield per hive fell by 50% between 2012 and 2022, whilst the number of occupied hives decreased by 44%. Beekeepers ranked climate change and declining flower abundance as the most important drivers of the decline. This raises concern for the future food and economic security of this region, where honey sales contribute to 16% of total household income, and where Apis cerana cerana plays a major role in crop pollination, contributing more than 50% of all flower visits to apple, cucumber, and pumpkin. To mitigate further declines, we promote native habitat and wildflower preservation, and using well-insulated log hives to buffer bees against the increasingly extreme temperature fluctuations.

Different currencies for calculating resource phenology result in opposite inferences about trophic mismatches.

Shifts in phenology are among the key responses of organisms to climate change. When rates of phenological change differ between interacting species they may result in phenological asynchrony. Studies have found conflicting patterns concerning the direction and magnitude of changes in synchrony, which have been attributed to biological factors. A hitherto overlooked additional explanation are differences in the currency used to quantify resource phenology, such as abundance and biomass. Studying an insectivorous bird (the sanderling) and its prey, we show that the median date of cumulative arthropod biomass occurred, on average, 6.9 days after the median date of cumulative arthropod abundance. In some years this difference could be as large as 21 days. For 23 years, hatch dates of sanderlings became less synchronized with the median date of arthropod abundance, but more synchronized with the median date of arthropod biomass. The currency-specific trends can be explained by our finding that mean biomass per arthropod specimen increased with date. Using a conceptual simulation, we show that estimated rates of phenological change for abundance and biomass can differ depending on temporal shifts in the size distribution of resources. We conclude that studies of trophic mismatch based on different currencies for resource phenology can be incompatible with each other.

Optimizing insect metabarcoding using replicated mock communities.

1: Metabarcoding (high-throughput sequencing of marker gene amplicons) has emerged as a promising and cost-effective method for characterizing insect community samples. Yet, the methodology varies greatly among studies and its performance has not been systematically evaluated to date. In particular, it is unclear how accurately metabarcoding can resolve species communities in terms of presence-absence, abundances, and biomass. 2: Here we use mock community experiments and a simple probabilistic model to evaluate the effect of different DNA extraction protocols on metabarcoding performance. Specifically, we ask four questions: (Q1) How consistent are the recovered community profiles across replicate mock communities?; (Q2) How does the choice of lysis buffer affect the recovery of the original community?; (Q3) How are community estimates affected by differing lysis times and homogenization?; and (Q4) Is it possible to obtain adequate species abundance estimates through the use of biological spike-ins? 3: We show that estimates are quite variable across community replicates. In general, a mild lysis protocol is better at reconstructing species lists and approximate counts, while homogenization is better at retrieving biomass composition. Small insects are more likely to be detected in lysates, while some tough species require homogenization to be detected. Results are less consistent across biological replicates for lysates than for homogenates. Some species are associated with strong PCR amplification bias, which complicates the reconstruction of species counts. Yet, with adequate spike-in data, species abundance can be determined with roughly 40% standard error for homogenates, and with roughly 50% standard error for lysates, under ideal conditions. In the latter case, however, this often requires species-specific reference data, while spike-in data generalizes better across species for homogenates. 4: We conclude that a non-destructive, mild lysis approach shows the highest promise for presence/absence description of the community, while also allowing future morphological or molecular work on the material. However, homogenization protocols perform better for characterizing community composition, in particular in terms of biomass.

The role of seasonality in shaping the interactions of honeybees with other taxa.

The Eltonian niche of a species is defined as its set of interactions with other taxa. How this set varies with biotic, abiotic and human influences is a core question of modern ecology. In seasonal environments, the realized Eltonian niche is likely to vary due to periodic changes in the occurrence and abundance of interaction partners and changes in species behavior and preferences. Also, human management decisions may leave strong imprints on species interactions. To compare the impact of seasonality to that of management effects, honeybees provide an excellent model system. Based on DNA traces of interaction partners archived in honey, we can infer honeybee interactions with floral resources and microbes in the surrounding habitats, their hives, and themselves. Here, we resolved seasonal and management-based impacts on honeybee interactions by sampling beehives repeatedly during the honey-storing period of honeybees in Finland. We then use a genome-skimming approach to identify the taxonomic contents of the DNA in the samples. To compare the effects of the season to the effects of location, management, and the colony itself in shaping honeybee interactions, we used joint species distribution modeling. We found that honeybee interactions with other taxa varied greatly among taxonomic and functional groups. Against a backdrop of wide variation in the interactions documented in the DNA content of honey from bees from different hives, regions, and beekeepers, the imprint of the season remained relatively small. Overall, a honey-based approach offers unique insights into seasonal variation in the identity and abundance of interaction partners among honeybees. During the summer, the availability and use of different interaction partners changed substantially, but hive- and taxon-specific patterns were largely idiosyncratic as modified by hive management. Thus, the beekeeper and colony identity are as important determinants of the honeybee's realized Eltonian niche as is seasonality.

Site-specific length-biomass relationships of arctic arthropod families are critical for accurate ecological inferences.

Arthropods play a crucial role in terrestrial ecosystems, for instance in mediating energy fluxes and in forming the food base for many organisms. To better understand their functional role in such ecosystem processes, monitoring of trends in arthropod biomass is essential. Obtaining direct measurements of the body mass of individual specimens is laborious. Therefore, these data are often indirectly acquired by utilizing allometric length-biomass relationships based on a correlative parameter, such as body length. Previous studies have often used such relationships with a low taxonomic resolution and/or small sample size and/or adopted regressions calibrated in different biomes. Despite the scientific interest in the ecology of arctic arthropods, no site-specific family-level length-biomass relationships have hitherto been published. Here we present 27 family-specific length-biomass relationships from two sites in the High Arctic: Zackenberg in northeast Greenland and Knipovich in north Taimyr, Russia. We show that length-biomass regressions from different sites within the same biome did not affect estimates of phenology but did result in substantially different estimates of arthropod biomass. Estimates of daily biomass at Zackenberg were on average 24% higher when calculated using regressions for Knipovich compared to using regressions for Zackenberg. In addition, calculations of daily arthropod biomass at Zackenberg based on order-level regressions from frequently cited studies in literature revealed overestimations of arthropod biomass ranging from 69.7% to 130% compared to estimates based on regressions for Zackenberg. Our results illustrate that the use of allometric relationships from different sites can significantly alter the biological interpretation of, for instance, the interaction between insectivorous birds and their arthropod prey. We conclude that length-biomass relationships should be locally established rather than being based on global relationships.

Mixed effects of a national protected area network on terrestrial and freshwater biodiversity.

Protected areas are considered fundamental to counter biodiversity loss. However, evidence for their effectiveness in averting local extinctions remains scarce and taxonomically biased. We employ a robust counterfactual multi-taxon approach to compare occupancy patterns of 638 species, including birds (150), mammals (23), plants (39) and phytoplankton (426) between protected and unprotected sites across four decades in Finland. We find mixed impacts of protected areas, with only a small proportion of species explicitly benefiting from protection-mainly through slower rates of decline inside protected areas. The benefits of protection are enhanced for larger protected areas and are traceable to when the sites were protected, but are mostly unrelated to species conservation status or traits (size, climatic niche and threat status). Our results suggest that the current protected area network can partly contribute to slow down declines in occupancy rates, but alone will not suffice to halt the biodiversity crisis. Efforts aimed at improving coverage, connectivity and management will be key to enhance the effectiveness of protected areas towards bending the curve of biodiversity loss.

Diverging trends and drivers of Arctic flower production in Greenland over space and time.

The Arctic is warming at an alarming rate. While changes in plant community composition and phenology have been extensively reported, the effects of climate change on reproduction remain poorly understood. We quantified multidecadal changes in flower density for nine tundra plant species at a low- and a high-Arctic site in Greenland. We found substantial changes in flower density over time, but the temporal trends and drivers of flower density differed both between species and sites. Total flower density increased over time at the low-Arctic site, whereas the high-Arctic site showed no directional change. Within and between sites, the direction and rate of change differed among species, with varying effects of summer temperature, the temperature of the previous autumn and the timing of snowmelt. Finally, all species showed a strong trade-off in flower densities between successive years, suggesting an effective cost of reproduction. Overall, our results reveal region- and taxon-specific variation in the sensitivity and responses of co-occurring species to shared climatic drivers, and a clear cost of reproductive investment among Arctic plants. The ultimate effects of further changes in climate may thus be decoupled between species and across space, with critical knock-on effects on plant species dynamics, food web structure and overall ecosystem functioning. Supplementary Information: The online version contains supplementary material available at 10.1007/s00300-023-03164-2.

Little directional change in the timing of Arctic spring phenology over the past 25 years.

With the global change in climate, the Arctic has been pinpointed as the region experiencing the fastest rates of change. As a result, Arctic biological responses-such as shifts in phenology-are expected to outpace those at lower latitudes. 15 years ago, a decade-long dataset from Zackenberg in High Arctic Greenland revealed rapid rates of phenological change.1 To explore how the timing of spring phenology has developed since, we revisit the Zackenberg time series on flowering plants, arthropods, and birds. Drawing on the full 25-year period of 1996-2020, we find little directional change in the timing of events despite ongoing climatic change. We attribute this finding to a shift in the temporal patterns of climate conditions, from previous directional change to current high inter-annual variability. Additionally, some taxa appear to have reached the limits of their phenological responses, resulting in a leveling off in their phenological responses in warm years. Our findings demonstrate the importance of long-term monitoring of taxa from across trophic levels within the community, allowing for detecting shifts in sensitivities and responses and thus for updated inference in the light of added information.

FAVIS: Fast and versatile protocol for non-destructive metabarcoding of bulk insect samples.

Insects are diverse and sustain essential ecosystem functions, yet remain understudied. Recent reports about declines in insect abundance and diversity have highlighted a pressing need for comprehensive large-scale monitoring. Metabarcoding (high-throughput bulk sequencing of marker gene amplicons) offers a cost-effective and relatively fast method for characterizing insect community samples. However, the methodology applied varies greatly among studies, thus complicating the design of large-scale and repeatable monitoring schemes. Here we describe a non-destructive metabarcoding protocol that is optimized for high-throughput processing of Malaise trap samples and other bulk insect samples. The protocol details the process from obtaining bulk samples up to submitting libraries for sequencing. It is divided into four sections: 1) Laboratory workspace preparation; 2) Sample processing-decanting ethanol, measuring the wet-weight biomass and the concentration of the preservative ethanol, performing non-destructive lysis and preserving the insect material for future work; 3) DNA extraction and purification; and 4) Library preparation and sequencing. The protocol relies on readily available reagents and materials. For steps that require expensive infrastructure, such as the DNA purification robots, we suggest alternative low-cost solutions. The use of this protocol yields a comprehensive assessment of the number of species present in a given sample, their relative read abundances and the overall insect biomass. To date, we have successfully applied the protocol to more than 7000 Malaise trap samples obtained from Sweden and Madagascar. We demonstrate the data yield from the protocol using a small subset of these samples.

Convergence of dominance and neglect in flying insect diversity.

Most of arthropod biodiversity is unknown to science. Consequently, it has been unclear whether insect communities around the world are dominated by the same or different taxa. This question can be answered through standardized sampling of biodiversity followed by estimation of species diversity and community composition with DNA barcodes. Here this approach is applied to flying insects sampled by 39 Malaise traps placed in five biogeographic regions, eight countries and numerous habitats (>225,000 specimens belonging to >25,000 species in 458 families). We find that 20 insect families (10 belonging to Diptera) account for >50% of local species diversity regardless of clade age, continent, climatic region and habitat type. Consistent differences in family-level dominance explain two-thirds of variation in community composition despite massive levels of species turnover, with most species (>97%) in the top 20 families encountered at a single site only. Alarmingly, the same families that dominate insect diversity are 'dark taxa' in that they suffer from extreme taxonomic neglect, with little signs of increasing activities in recent years. Taxonomic neglect tends to increase with diversity and decrease with body size. Identifying and tackling the diversity of 'dark taxa' with scalable techniques emerge as urgent priorities in biodiversity science.

Food webs coupled in space: Consumer foraging movement affects both stocks and fluxes.

The exchange of material and individuals between neighboring food webs is ubiquitous and affects ecosystem functioning. Here, we explore animal foraging movement between adjacent, heterogeneous habitats and its effect on a suite of interconnected ecosystem functions. Combining dynamic food web models with nutrient-recycling models, we study foraging across habitats that differ in fertility and plant diversity. We found that net foraging movement flowed from high to low fertility or high to low diversity and boosted stocks and flows across the whole loop of ecosystem functions, including biomass, detritus, and nutrients, in the recipient habitat. Contrary to common assumptions, however, the largest flows were often between the highest and intermediate fertility habitats rather than highest and lowest. The effect of consumer influx on ecosystem functions was similar to the effect of increasing fertility. Unlike fertility, however, consumer influx caused a shift toward highly predator-dominated biomass distributions, especially in habitats that were unable to support predators in the absence of consumer foraging. This shift resulted from both direct and indirect effects propagated through the interconnected ecosystem functions. Only by considering both stocks and fluxes across the whole loop of ecosystem functions do we uncover the mechanisms driving our results. In conclusion, the outcome of animal foraging movements will differ from that of dispersal and diffusion. Together we show how considering active types of animal movement and the interconnectedness of ecosystem functions can aid our understanding of the patchy landscapes of the Anthropocene.

Different spatial structure of plant-associated fungal communities above- and belowground.

The distribution and community assembly of above- and belowground microbial communities associated with individual plants remain poorly understood, despite its consequences for plant-microbe interactions and plant health. Depending on how microbial communities are structured, we can expect different effects of the microbial community on the health of individual plants and on ecosystem processes. Importantly, the relative role of different factors will likely differ with the scale examined. Here, we address the driving factors at a landscape level, where each individual unit (oak trees) is accessible to a joint species pool. This allowed to quantify the relative effect of environmental factors and dispersal on the distribution of two types of fungal communities: those associated with the leaves and those associated with the soil of Quercus robur trees in a landscape in southwestern Finland. Within each community type, we compared the role of microclimatic, phenological, and spatial variables, and across community types, we examined the degree of association between the respective communities. Most of the variation in the foliar fungal community was found within trees, whereas soil fungal community composition showed positive spatial autocorrelation up to 50 m. Microclimate, tree phenology, and tree spatial connectivity explained little variation in the foliar and soil fungal communities. Foliar and soil fungal communities differed strongly in community structure, with no significant concordance detected between them. We provide evidence that foliar and soil fungal communities assemble independent of each other and are structured by different ecological processes.

Stable pollination service in a generalist high Arctic community despite the warming climate.

Insects provide key pollination services in most terrestrial biomes, but this service depends on a multistep interaction between insect and plant. An insect needs to visit a flower, receive pollen from the anthers, move to another conspecific flower, and finally deposit the pollen on a receptive stigma. Each of these steps may be affected by climate change, and focusing on only one of them (e.g., flower visitation) may miss important signals of change in service provision. In this study, we combine data on visitation, pollen transport, and single-visit pollen deposition to estimate functional outcomes in the high Arctic plant-pollinator network of Zackenberg, Northeast Greenland, a model system for global warming-associated impacts in pollination services. Over two decades of rapid climate warming, we sampled the network repeatedly: in 1996, 1997, 2010, 2011, and 2016. Although the flowering plant and insect communities and their interactions varied substantially between years, as expected based on highly variable Arctic weather, there was no detectable directional change in either the structure of flower-visitor networks or estimated pollen deposition. For flower-visitor networks compiled over a single week, species phenologies caused major within-year variation in network structure despite consistency across years. Weekly networks for the middle of the flowering season emerged as especially important because most pollination service can be expected to be provided by these large, highly nested networks. Our findings suggest that pollination ecosystem service in the high Arctic is remarkably resilient. This resilience may reflect the plasticity of Arctic biota as an adaptation to extreme and unpredictable weather. However, most pollination service was contributed by relatively few fly taxa (Diptera: Spilogona sanctipauli and Drymeia segnis [Muscidae] and species of Rhamphomyia [Empididae]). If these key pollinators are negatively affected by climate change, network structure and the pollination service that depends on it would be seriously compromised.

Emerging technologies revolutionise insect ecology and monitoring.

Insects are the most diverse group of animals on Earth, but their small size and high diversity have always made them challenging to study. Recent technological advances have the potential to revolutionise insect ecology and monitoring. We describe the state of the art of four technologies (computer vision, acoustic monitoring, radar, and molecular methods), and assess their advantages, current limitations, and future potential. We discuss how these technologies can adhere to modern standards of data curation and transparency, their implications for citizen science, and their potential for integration among different monitoring programmes and technologies. We argue that they provide unprecedented possibilities for insect ecology and monitoring, but it will be important to foster international standards via collaboration.

Reconstructing the ecosystem context of a species: Honey-borne DNA reveals the roles of the honeybee.

To assess a species' impact on its environment-and the environment's impact upon a species-we need to pinpoint its links to surrounding taxa. The honeybee (Apis mellifera) provides a promising model system for such an exercise. While pollination is an important ecosystem service, recent studies suggest that honeybees can also provide disservices. Developing a comprehensive understanding of the full suite of services and disservices that honeybees provide is a key priority for such a ubiquitous species. In this perspective paper, we propose that the DNA contents of honey can be used to establish the honeybee's functional niche, as reflected by ecosystem services and disservices. Drawing upon previously published genomic data, we analysed the DNA found within 43 honey samples from Northern Europe. Based on metagenomic analysis, we find that the taxonomic composition of DNA is dominated by a low pathogenicity bee virus with 40.2% of the reads, followed by bacteria (16.7%), plants (9.4%) and only 1.1% from fungi. In terms of ecological roles of taxa associated with the bees or taxa in their environment, bee gut microbes dominate the honey DNA, with plants as the second most abundant group. A range of pathogens associated with plants, bees and other animals occur frequently, but with lower relative read abundance, across the samples. The associations found here reflect a versatile the honeybee's role in the North-European ecosystem. Feeding on nectar and pollen, the honeybee interacts with plants-in particular with cultivated crops. In doing so, the honeybee appears to disperse common pathogens of plants, pollinators and other animals, but also microbes potentially protective of these pathogens. Thus, honey-borne DNA helps us define the honeybee's functional niche, offering directions to expound the benefits and drawbacks of the associations to the honeybee itself and its interacting organisms.

Beyond body size-new traits for new heights in trait-based modelling of predator-prey dynamics.

Food webs map feeding interactions among species, providing a valuable tool for understanding and predicting community dynamics. Using species' body sizes is a promising avenue for parameterizing food-web models, but such approaches have not yet been able to fully recover observed community dynamics. Such discrepancies suggest that traits other than body size also play important roles. For example, differences in species' use of microhabitat or non-consumptive effects of intraguild predators may affect dynamics in ways not captured by body size. In Laubmeier et al. (2018), we developed a dynamic food-web model incorporating microhabitat and non-consumptive predator effects in addition to body size, and used simulations to suggest an optimal sampling design of a mesocosm experiment to test the model. Here, we perform the mesocosm experiment to generate empirical time-series of insect herbivore and predator abundance dynamics. We minimize least squares error between the model and time-series to determine parameter values of four alternative models, which differ in terms of including vs excluding microhabitat use and non-consumptive predator-predator effects. We use both statistical and expert-knowledge criteria to compare the models and find including both microhabitat use and non-consumptive predator-predator effects best explains observed aphid and predator population dynamics, followed by the model including microhabitat alone. This ranking suggests that microhabitat plays a larger role in driving population dynamics than non-consumptive predator-predator effects, although both are clearly important. Our results illustrate the importance of additional traits alongside body size in driving trophic interactions. They also point to the need to consider trophic interactions and population dynamics in a wider community context, where non-trophic impacts can dramatically modify the interplay between multiple predators and prey. Overall, we demonstrate the potential for utilizing traits beyond body size to improve trait-based models and the value of iterative cycling between theory, data and experiment to hone current insights into how traits affect food-web dynamics.

Comparison of traditional and DNA metabarcoding samples for monitoring tropical soil arthropods (Formicidae, Collembola and Isoptera).

The soil fauna of the tropics remains one of the least known components of the biosphere. Long-term monitoring of this fauna is hampered by the lack of taxonomic expertise and funding. These obstacles may potentially be lifted with DNA metabarcoding. To validate this approach, we studied the ants, springtails and termites of 100 paired soil samples from Barro Colorado Island, Panama. The fauna was extracted with Berlese-Tullgren funnels and then either sorted with traditional taxonomy and known, individual DNA barcodes ("traditional samples") or processed with metabarcoding ("metabarcoding samples"). We detected 49 ant, 37 springtail and 34 termite species with 3.46 million reads of the COI gene, at a mean sequence length of 233 bp. Traditional identification yielded 80, 111 and 15 species of ants, springtails and termites, respectively; 98%, 37% and 100% of these species had a Barcode Index Number (BIN) allowing for direct comparison with metabarcoding. Ants were best surveyed through traditional methods, termites were better detected by metabarcoding, and springtails were equally well detected by both techniques. Species richness was underestimated, and faunal composition was different in metabarcoding samples, mostly because 37% of ant species were not detected. The prevalence of species in metabarcoding samples increased with their abundance in traditional samples, and seasonal shifts in species prevalence and faunal composition were similar between traditional and metabarcoding samples. Probable false positive and negative species records were reasonably low (13-18% of common species). We conclude that metabarcoding of samples extracted with Berlese-Tullgren funnels appear suitable for the long-term monitoring of termites and springtails in tropical rainforests. For ants, metabarcoding schemes should be complemented by additional samples of alates from Malaise or light traps.

Elevation and plant species identity jointly shape a diverse arbuscular mycorrhizal fungal community in the High Arctic.

Knowledge about the distribution and local diversity patterns of arbuscular mycorrhizal (AM) fungi are limited for extreme environments such as the Arctic, where most studies have focused on spore morphology or root colonization. We here studied the joint effects of plant species identity and elevation on AM fungal distribution and diversity. We sampled roots of 19 plant species in 18 locations in Northeast Greenland, using next generation sequencing to identify AM fungi. We studied the joint effect of plant species, elevation and selected abiotic conditions on AM fungal presence, richness and composition. We identified 29 AM fungal virtual taxa (VT), of which six represent putatively new VT. Arbuscular mycorrhizal fungal presence increased with elevation, and as vegetation cover and the active soil layer decreased. Arbuscular mycorrhizal fungal composition was shaped jointly by elevation and plant species identity. We demonstrate that the Arctic harbours a relatively species-rich and nonrandomly distributed diversity of AM fungi. Given the high diversity and general lack of knowledge exposed herein, we encourage further research into the diversity, drivers and functional role of AM fungi in the Arctic. Such insight is urgently needed for an area with some of the globally highest rates of climate change.